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A method is presented for solving the characteristic ini- 
tial value problem for the collision and subsequent nonlinear 
interaction of plane gravitational or gravitational and elec- 
tromagnetic waves in a Minkowski background. This method 
generalizes the monodromy transform approach to fields with 
nonanalytic behaviour on the characteristics inherent to waves 
with distinct wave fronts. The crux of the method is in a refor- 
mulation of the main nonlinear symmetry reduced field equa- 
tions as linear integral equations whose solutions are deter- 
mined by generalized ("dynamical") monodromy data which 
evolve from data specified on the initial characteristics (the 
wavefronts) . 
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I. INTRODUCTION 

The collision and subsequent nonlinear interaction be- 
tween plane gravitational or gravitational and electro- 
magnetic waves propagating with distinct wavefronts in 
a Minkowski background is a well formulated character- 
istic initial value problem. However, even the discovery 
of the integrability of the main field equations for this 
situation did not lead to a solution of this complex non- 
linear problem. None of the existing solution generating 
methods have been found suitable for an effective con- 
struction of the solution in the wave interaction region 
starting from given characteristic initial data determined 
by the approaching waves. 

The structure of the governing field equations for col- 
liding plane waves, physical and geometrical interpreta- 
tions, and various particular solutions and techniques, 
have been described in [fij. For colliding plane gravita- 
tional waves with aligned constant polarizations, the vac- 
uum Einstein equations are reducible to the linear Eulcr- 
Poisson-Darboux equation. In this case, the correspond- 
ing characteristic initial value problem can be solved us- 
ing the generalized version of Abel's transform ||. 

However, when the polarizations of the approaching 
gravitational waves are not constant and aligned, or in 
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the presence of electromagnetic waves, the governing 
equations are essentially nonlinear and are equivalent to 
the hyperbolic form of the Ernst equations. For this case, 
an appreciable number of particular solutions are known. 
These have been found using the "inverse" method in 
which a formal solution in the interaction region is first 
constructed and the corresponding characteristic initial 
data for the approaching waves is only determined sub- 
sequently. Recently, infinite hierarchies of exact vacuum 
and electrovacuum solutions with an arbitrary number 
of free parameters were found || , and many of these are 
of the type appropriate for colliding plane waves. How- 
ever, it is not a simple technical task to simplify these 
solutions for particular cases and to calculate the corre- 
sponding characteristic initial data. 

For the analysis of the characteristic initial value prob- 
lem for the vacuum hyperbolic Ernst equation, Hauser 
and Ernst Q] have generalized their group-theoretical 
approach (which had been developed earlier for station- 
ary axisymmctric fields) and constructed a homogeneous 
Hilbert problem with corresponding matrix linear inte- 
gral equations. Many aspects of this problem, including 
the existence and uniqueness of solutions and a detailed 
proof of the Geroch conjecture, were elaborated in ||. 

However, a general scheme for the solution of various 
nonlinear initial and boundary value problems for inte- 
grable reductions of Einstein's equations had been de- 
veloped in the framework of the monodromy transform 
approach ||- ||. This can be applied to both the char- 
acteristic initial value problem and the Cauchy problem 
for the hyperbolic case, as well as some boundary prob- 
lems for the elliptic case. In this approach, every solution 
can be characterized by a set of functions of an auxiliary 
(spectral) parameter. These functions are interpreted as 
the monodromy data on the spectral plane of the fun- 
damental solution of an auxiliary overdetermined linear 
system associated with the (nonlinear) field equations. 
These monodromy data are nonevolving (i.e. coordinate 
independent) and, generally, can be chosen arbitrarily or 
specified in accordance with the properties of the solu- 
tion being sought. In particular, these data can be deter- 
mined (at least in principle) from the initial or bound- 
ary data. In this scheme, the solution of the initial or 
boundary value problem is determined by the solution 
of some linear singular integral equations whose kernel is 
constructed using these (specified) monodromy data. 

Recently a further generalization of this approach was 



1 



derived in |L0|], in which a new linear integral equation 
form of various hyperbolic integrable reductions of Ein- 
stein equations was constructed. The scalar kernel of 
these quasi-Fredholm equations depends on the mon- 
odromy data in a different way - referred to as "dynami- 
cal monodromy data" . Unlike the previous method, these 
data evolve, and their evolution is prescribed by the char- 
acteristic initial conditions. It appears that these new 
integral equations are better adapted to the construction 
of an effective solution of initial value problems because 
their coefficients carry more explicit information about 
the characteristic initial data and the corresponding an- 
alytical structures of the solutions desired. 

In this paper, we present a new approach to the so- 
lution of the colliding plane wave problem and also de- 
scribe a method that can be implemented in practice to 
derive explicit solutions. Our construction generalizes 
the monodromy transform approach whose formulation 
in the above mentioned papers was applicable only for 
fields that are analytically dependent on some special set 
of geometrically defined space-time coordinates. (One 
of these coordinates determines the measure of an area 
on the two-dimensional orbits of the space-time isome- 
try group. The other is its harmonic conjugate.) How- 
ever, for colliding plane waves propagating with distinct 
wave fronts on a Minkowski background, this analyticity 
is obviously violated on the wave fronts. Moreover, the 
physically accepted matching conditions (the "colliding 
plane wave conditions") imply regular behaviour of the 
field components near the wave fronts in some appropri- 
ate null coordinates [jlj . In terms of geometrically defined 
coordinates, this regularity leads to a specific singular 
behaviour of the coefficients of the associated linear sys- 
tem on these hypersurfaces, where the first derivatives of 
the field components become infinite. This gives rise to 
crucial consequences for the previous formulation. Addi- 
tional singularities appear on the spectral plane for some 
auxiliary functions. Also, the integral representations 
diverge at the singular points, while the normalization 
of the solutions remains a necessity. Another important 
phenomenon which arises in these singular cases is that 
the nonevolving monodromy data, which continue to ex- 
ist, lose their the most important property - their unam- 
biguous characterization of the solution. 

The solution of these problems arises from our recent 
observation that one of two linear integral "evolution 
equation" forms of the field equations derived in |fCf| ad- 
mits a generalization to the singular case. It is also im- 
portant that the dynamical monodromy data can still be 
used in this case to characterize the solutions. In this 
paper we present the generalized linear quasi-Fredholm 
integral "evolution equation" , which covers the singular 
case and opens a direct way for the construction of solu- 
tions for colliding plane waves from given characteristic 
initial data. We also briefly discuss some applications. 



II. ASSOCIATED LINEAR SYSTEM WITH 
SPECTRAL PARAMETER 

We base our construction on the Kinnersley-like 
overdetermined linear system with a free complex pa- 
rameter w for aJVxJV matrix function \I/(£, rj, w) whose 
integrability conditions are equivalent to the hyperbolic 
space-time symmetry reduction for N = 2 of the vacuum 
Einstein equations and for N = 3 of the electrovacuum 
Einstein-Maxwell equations - 

f 2i(w-£)<%* = U(£,r7)*, rankU=l, trU = i, 
\ 2i(w - rj)d v ^ = V(£,77)*, rank V = 1, trV = i 

in which the null coordinates £, r\ are certain linear com- 
binations of geometrically defined nonnull coordinates 
mentioned in the introduction. These equations must 
be supplemented with the additional constraints on its 
matrix integral |t]|| 

fftW* = W H <9W 

\ Wj( w ) = Wo (to), ~ 1 

where t denotes Hermitian conjugation (Wj(w) = 

Wq(w)) and tt is a constant matrix. For N = 2,3, 
the only nonzero components of Q are ft 12 = 1 and 
n 21 = -1. For N = 3, the condition W 33 = I should 
be also satisfied. The real null coordinates £, rj are cer- 
tain linear combinations of geometrically defined nonnull 
coordinates, say a, /3, mentioned in the introduction. 

For any solution of these conditions, the components 
of U, V and W can be identified with certain metric 
components and the electromagnetic potential and their 
derivatives. Without loss of generality, we impose the 
normalization conditions at the point denoted by (£ x , rj x ) 
at which the waves collide. These are ^f(£, x ,ri x ,w) = I 
and W (w) = Ai(w— /3 x )fi+diag (4a x , 4, 1), where a x — 
(£x — r)y)/2 and f3 x = (£ x + r] x )/2. A shift of origin and 
a rescaling of the coordinates £ and rj allow us to specify 
below £ x = 1 and r) x = — 1, so that a x = 1, (3 X = 0. 

III. THE COLLIDING PLANE WAVE PROBLEM 

For plane gravitational or gravitational and electro- 
magnetic waves with distinct wavefronts which collide 
in a Minkowski background, it is well known Jlj that 
there exist global null coordinates (it, v) such that the 
O'Brien-Singe matching conditions imposed on the wave- 
fronts u — and v — lead to a well posed characteristic 
initial value problem. However, the regularity of this 
problem in these global coordinates in the interaction re- 
gion (u > 0, v > 0) implies special relations between 
these coordinates and the coordinates (£, rj) which are 
nonanalytical on the boundaries. It is possible to use the 
coordinate freedom to put 

£ = I-2u n +, r) = -i + 2v n -, 
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where n± > 2 are specified as part of the initial data. 
Thus, the important specific of this characteristic initial 
value problem is that, in terms of geometrically defined 
coordinates (£,??), the first derivatives of the field com- 
ponents should be discontinuous and even unbounded on 
the wavefronts and at the point of collision. 

IV. THE ANALYTICAL STRUCTURE OF * ON 
THE SPECTRAL PLANE 

Everywhere below, \& denotes the fundamental so- 
lution of the associated linear system described above 
which is normalized at the point of collision [u = O^v = 0, 
i.e. (£ = l,r] = —1). As in the regular case (7]|J, the 
structure of the coefficients of the associated linear sys- 
tem for implies that the normalized fundamental solu- 
tion \&(£, rj, w) and its inverse possess four branch points 
on the spectral plane w, namely at w = £, w = rj, w = 1 
and to = — 1. The order of these points and our choice of 
the cuts L± joining them are indicated in the figure 

L_ L + 



Near these singular points, and on the cuts L±, the com- 
ponents of \I/ and its inverse (as functions of w for given 
u > and v> 0) possess in general the local structure: 

9(t,r),w) =$±(£,T],w) ®k±M + M±(f,»7,io) 

where, as in the regular case, the coordinate indepen- 
dent components of the row vectors k±(u>) constitute 
the "projective vectors" of the monodromy data. Their 
components, as well as the components of the matrices 
M±(£, 77, to), are regular (holomorphic) on the cuts L± 
with the corresponding index. In the analytical case 

we always have $ + (£,r),w) = yfijjjE^ty +(£> V, w ) and 

t|> - j Vi w ) = *J ^r^4' - (£j 7 h w ) where the vectors ip + 
and \p _ are holomorphic on the cuts L + and L_ respec- 
tively. In^the general case, the components of the column 
vectors ip ±(£,77,10) also have branch points at the end- 
points of the corresponding cuts L±, but the character 
of their singularities at the points 10 = 1 and w = — 1 
respectively (as well as the powers n + and 71 _ defined 
above) are determined by the initial data. 

V. THE INTEGRAL "EVOLUTION" EQUATIONS 
AND THE SOLUTION OF THE PROBLEM 

To set up the characteristic initial value problem, we 
introduce two matrix functions which are the normalized 
fundamental solutions of the linear ordinary differential 



equations - the restrictions of the associated linear sys- 
tem to the characteristics £ = 1 and 77 = — 1: 

2i(w - £)%*+ = U(£, -1) • *+, to) = I 

2i(w - 7])d v ^ - = V(l,ry) • *_(-l,ui) = I 

in which the coefficients are determined by the initial 
data for the fields on the corresponding characteristics. 
These matrices should be the characteristic initial data 
for the required solution for ^(£,77,10): 

*+(£,to) = #(£,-l,to), 9-(ji,w) = *(l,7?,to) 

The analytical structures of \t± on the spectral plane 
are very similar to those of \E'. Namely, *&±(w = 00) = I, 
is holomorphic outside L + and M>_ outside L_, and 
their local structures on these cuts are given by 

L+ : *+(£, w) = ij> 0+ (£, w) (g> k+(to) + M 0+ (£, w), 
L_ : *-(?7, w) = ip _(?7, 10) ® k_(w) + M -(?7, w) 

where k±(w) are the same as for and Mo+(£,to) and 
Mo_ (77, to) are holomorphic on L + and L_ respectively. 

We now introduce, in analogy with the regular case 
]l0| , the "evolution" or "scattering" matrices x±(£,?7,w), 
representing \&(£, 77, to) in two alternative forms 

¥(£,ri,w) = X + (t,ri,w) ■ *+(£,w) 
®(€,V,w) = X-(^V,w) ■ *_(7?,to). 

The crucial point is that the components of the ma- 
trix x + (£: 7 h w ) are holomorphic on L + and possess a 
jump on L_ only, while the components of the matrix 
X-(£> "Hi w) are holomorphic on L_ and possess a jump 
on L+. From the above it is clear also, that these jumps 
are represented by highly degenerate matrices and that 
X±(£> 7 7j ?w = 00 ) = I- These properties permit us to 
represent x± as Cauchy integrals of the form 

X + U,?7,w =I+— / r d£- 

7TZ J 10 - C+ 

£+ 

where [. . .] is the jump (a half of the difference between 
the left and right limits) of a function on a cut and 

m_(£,to)=k_(to)-*; 1 (£,tt;) 

m + (?7, to) = k+(io) • ^'Z 1 (r7, to) 

represent a new (evolving) kind of monodromy data in- 
troduced for the regular case in [|l0| and called there the 
"dynamical" monodromy data. It is necessary to note 
here our conjecture of the convergency of the Cauchy in- 
tegrals for x± at w = — 1 and to = 1 respectively which 
is confirmed for particular examples for any n± > 2. 
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The above alternative representations for St should sat- 
isfy an obvious condition X + *+ = X-*-- This condi- 
tion considered on the cuts L± together with the con- 
structed integral representations for x ± leads to the lin- 
ear integral equations (for brevity we omit here the para- 
metric dependence of all objects upon £ and 77): 

<t>+(r+) - J S+(T + ,C-)<S>-((-)d(- = * +(r+) 
z,_ 



*-(r_)-y S_(t_,C+)cMC + K + = c|>o-(t-) 

where t + ,£ + 6 L + , t_,£_ G L_ and the vector func- 
tions cjj +(t + ), c)j -(t_) and their initial values cf> o+(t + ), 
cj)o-(r-) are the jumps [tp+] T+ , [tp _]r_ and [^ +]r+, 
[4* o-]r_ respectively. The scalar kernels are given by 
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i7r(C- - r+) 

S-(n,T-,C+) = —77 7(m + (?7,C+) • $>o-(ri,T- 

Z7T(C+ -t_) \ 

As in the regular case, the coefficients of the above inte- 
gral equations are determined by the initial data. How- 
ever, these generalized "evolution equations" possess a 
more complicated singular structure. They also can be 
easily decoupled into two independent equations for cf> + 
and cf) _, but we omit these calculations here. 

Our construction of solutions for colliding plane waves 
begins with the constants n± > 2 which determine the 
degree of non-smoothness of the fields on the wavefronts 
and the characteristic initial data for the fields in terms 
of the Ernst potentials £(u, v), &(u, v) which characterise 
every solution. These data, viz. £ + {u), <&+(u), £-(«), 
and <£>_(u), should be chosen to satisfy the normalization 
conditions £±(0) = —1, $±(0) = and two wavefront 
regularity conditions 0: |£±(0)| 2 + 4|$' ± (0)| 2 = 8(1 - 
l/n±). With this data we have to solve the ordinary 
differential equations for *f?± and the integral equations 
for cf) ± . The Ernst potentials of the sought-for solution 
can be evaluated then as 

£{u,v)=£+{u) + l J ( ei -(Jj_(C-))(m_(C-)-e 2 K_ 
= £_{v) + lj (e 1 .cMC+))(m + (C+)-e 2 )dC+ 

L+ 

$(u,«) =*+(«)-£ J (ei.*_(C-))(m_(C_)-e3)dC- 
= #_(«)- § J (ei • * + (C+))(m + (C+) • e 3 )dC+ 

where ei = {1, 0, 0}, e 2 = {0, 1, 0} and e 3 = {0, 0, 1}. 
As a simple test, we consider ri + = ?i_ = 2 and chose 

£+ = -(l-uf, £- = -(l-v) 2 



For these vacuum data our calculations lead immediately 
to the Khan-Penrose solution for the collision of impul- 
sive gravitational waves with colinear polarizations. Sim- 
ilar calculations for the initial data (5± are real constants) 



£+ = -l + 2e 



£- = -l + 2e l5 



give the Nutku-Halil solution for the collision of these 
waves with noncolinear polarizations. The initial data 



£+ = -1, $+=u, 



£- = -1, = v 



lead to the Bell-Szekeres solution for the collision of elec- 
tromagnetic step waves with aligned polarizations. When 
the polarizations of these waves is nonaligned, we change 



the data above to 
produces 



ve 
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, (7 is a real constant) . This 



£(u,v) = —l, <f>(it, v) — u\J 1 - v 2 + v\J 1 - u V 7 

which leads to a nontrivial (nondiagonal) metric (see 
]Tl[). Clearly, various new solutions can be found in the 
same way for different choices of the initial data. 
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